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Abstract. We study the characteristic exponents of the Hamil- 
tonian system of n (> 2) point masses mi, . . . , m n freely falling 
in the vertical half line {q\ q > 0} under constant gravitation and 
colliding with each other and the solid floor q = elastically. 
This model was introduced and first studied by M. Wojtkowski. 
Hereby we prove his conjecture: All relevant characteristic (Lya- 
punov) exponents of the above dynamical system are nonzero, 
provided that mi > • • • > m n (i. e. the masses do not increase as 
we go up) and mi ^ m 2 . 

1. Introduction 



In his paper [W(1990)-I] M. Wojtkowski introduced the following Hamiltonian 
dynamical system with discontinuities: There is a vertical half line {q\ q > 0} given 
and n (> 2) point particles with masses mi > rri2 > • • • > m n > and positions 
< qi < q2 < • • • < q n are moving on this half line so that they are subjected to a 
constant gravitational acceleration a = — 1 (they fall down), they collide elastically 
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with each other, and the first (lowest) particle also collides elastically with the hard 
floor q = 0. We fix the total energy 



by taking H = 1. The arising Hamiltonian flow with collisions (M, {il) t \ t E M} , n) 
(/j is the Liouville measure) is the subject of this paper. 

Before formulating the result of this article, however, it is worth mentioning here 
three important facts: 

(1) Since the phase space M is compact, the Liouville measure \i is finite. 

(2) The phase points x E M for which the orbit {^(x)] tel} hits at least one 
singularity (i. e. a multiple collision) are contained in a countable union of proper, 
smooth submanifolds of M and, therefore, such points form a set of \x measure zero. 

(3) For /x-almost every phase point x E M the collision moments of the orbit 
{ip t (x) \ t E R} do not have any finite accumulation point, see the Appendix. 

In the paper [W(1990)-I] Wojtkowski formulated his main conjecture pertaining 
to the dynamical system (M, {ip t \ t E R} , fi): 

Wojtkowski's Conjecture. If m\ > > ■ ■ ■ > m n > and mi ^ m n , then 
all but one characteristic (Lyapunov) exponents of the flow (M, t E R} , fi) are 
nonzero. 

Remarks. 

1. The only exceptional exponent zero must correspond to the flow direction. 

2. The condition of nonincreasing masses (as above) is essential for establishing 
the invariance of the symplectic cone field — an important condition for obtaining 
nonzero characteristic exponents. As Wojtkowski pointed out in Proposition 4 of 
[W(1990)-I], if n = 2 and mi < m2, then there exists a linearly stable periodic 
orbit, thus dimming the chances of proving ergodicity. 



In pursuing the goal of proving this conjecture, Wojtkowski obtained the follow- 
ing results in [W(1990)-I]: 
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Proposition 1. For every e > there is a 5 > such that if mi > mi > • • • > 

m n > and < 5, then (M, tel},/i) has exactly one zero charac- 

m\ 

teristic exponent except possibly on a set of [i measure < e. 

Proposition 2. If there are exactly I groups of particles with equal masses, I > 2, 
containing k±, . . . ,ki particles respectively, the greatest common divisor of ki, . . . ,ki 
is one and mi > mi > • ■ ■ > m n > 0, then {i^ 1 } has exactly one zero characteristic 
exponent on a set of positive Liouville measure. 

Proposition 3. If n = 3 and mi > m^ > m%, then has exactly one zero 

characteristic exponent ^-almost everywhere. 

In the subsequent article [W(1990)-II] Wojtkowski replaced the linear potential 
U(q) = q of constant gravitation by a varying gravitational force with potential 
U(q) for which U'(q) > and U"(q) < 0. (The usual gravitational potential 
U(q) = ^— belongs to this category!) He proved there that in the falling ball 

q + c 

system with such a potential U(q) all relevant characteristic exponents are nonzero 
almost everywhere. 

The result of this paper is a slightly weakened version of Wojtkowski' s conjecture: 

Theorem. // mi > mi > ms > • • • > m n > 0, then fi-almost everywhere all but 
one characteristic exponents of the flow (M, \ t G K} , y) are nonzero. 

We are closing this brief introduction by mentioning that in his work [Ch(1993)] 
N. I. Chernov significantly relaxed a condition of the Liverani-Wojtkowski local 
ergodicity theorem for symplectomorphisms, [L-W(1995)]. (This theorem is a gen- 
eralization of the celebrated Theorem on Local Ergodicity for semi-dispersing bil- 
liards by Chernov and Sinai, [S-Ch(1987)].) The ominous condition is the "proper 
alignment" of the singularity manifolds, Condition D in Section 7 of [L-W(1995)]. 
This condition is easily seen to be violated by the system of falling balls (see Section 
14.C of [L-W(1995)]), but the relaxed condition 5' of Chernov's paper [Ch(1993)] 
is very likely to hold for the system of falling balls. Thus, thanks to Chernov's im- 
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provement in the local ergodicity theorem, the ergodicity of the falling ball system 
essentially got within our reach! 

For a more detailed introduction to this subject, and for a thoroughly assembled 
collection of references and historical remarks, the reader is kindly referred to the 
introduction of the paper [W(1990)-I]. 

Remark. It is easy to see that the studying of the falling point particles on the 
vertical half line {q\ q > 0} is not a restriction of generality as compared to the 
systems of 1-D balls (hard rods) of length 2r. Namely, the simple change in the 
kinetic data qi i— > qi — (2i — l)r, Vi i— > Vi, H \— > H —r Yl7=i(^~ l)mj (the change of 
the fixed level of energy) establishes an isomorphism between the hard rod system 
and the point particle model. 



2. Prerequisites 



The upcoming brief survey of our dynamical system and the related technical- 
ities will narrowly follow the approach of sections 1-3 of [W(1990)-I]. A thorough 
description of the falling ball model and detailed references can be found in that 
article. 

We consider the following Hamiltonian system with discontinuities: Given the 
vertical half line {q\ q > 0}, n (> 2) point particles (one dimensional "balls" labelled 
by 1,2, ... ,n) with positions < q\ < q2 < • • • < q n are moving in that half line 
so that they fall down under a constant gravitational acceleration a = — 1, they 
collide with each other elastically whenever they hit each other, and the first particle 
bounces back from the hard floor q = elastically when it hits the floor. Denote 
by Vi = <ji the velocity, by pi = rriiVi the momentum, and by hi = rriiqi + -rriivf = 

2 

m i<li + tt— the energy of the i-th ball. (The quantity rriiqi is the potential energy 
Zrrii 

arising from the gravitation.) The extended phase space (without fixing the energy) 
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of this mechanical system is then 

(2.1) N = {(q,p) G R n x R n \ < q x < ■ ■ ■ < q n } . 
The manifold (with boundary) N carries the usual symplectic form 

n n 

(2.2) u) = dqi A dpi = dhi A dvi. 



i=i i=i 



The Hamiltonian function is H(q,p) = $^i=i ^ an d the arising Hamiltonian flow 
{ip t \ t G M} is determined by the usual formalism 



(2.3) 



. Pi dH 

nii dpi 

dH 

Pi = -mi = 



dqi 

(i = 1, . . . , n) between collisions. 

A collision of type (z, i + 1) (i = 1, . . . , n — 1) occurs when qi = qi+i. Then the 
velocities and the momenta of the colliding particles get transformed according to 
the law of elastic collisions: 



(2.4) 



( v+ = 7iur + (1 - n/i)vr+i 
v? +1 = (l + 7iK r -7*^i 
vt = HPi + (! + li)Pi + i 
Pt+i = (! - li)Pi ~ liPi+n 



ffl^ — 777, 

where 7$ = — — , see also (3) of [W(1990)-Il. (Here the superscript — (+) 

rrii + m i+1 ' 

refers to the kinetic data measured right before (after) the considered collision.) 
At a floor collision q\ = we obviously have 



(2.5) 



Pi = -Pi 



We introduce the following notations for the several components of the boundary 
<9N of N: 



(2.6) 



ONf = {(q,p) G N| q { = q i+1 k ± (v t - v i+1 ) < 0} 
9N± = U^dNf 
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for i = 0, . . . , n — 1, where, by convention, qo = vq = 0. Then the collision map 
$i : <9N~ -> <9N+ (i = 0, ... ,n - 1), determined by (2.4)-(2.5) and by the con- 
dition that the positions do not change at collisions, turns out to be a symplectic 
diffeomorphism of <9N~ onto dN+, see Section 2 of [W(1990)-I]. (Here the (2n- 1)- 
dimensional manifold <9N naturally inherits the pseudo-symplectic structure to |~<9N, 
i. e. the restriction of uj to <9N.) 

Now we fix the total energy by taking H = H = 1, and consider the restriction 
of the Hamiltonian flow with collisions t G IR} to the energy hypersurface 



The corresponding boundary components of M are denoted by dNl i and <9M . 

The (invariant) Liouville measure v in the extended phase space N is defined via 
the volume element V = cuA- • -Ao; (the n-th exterior power of ui). The corresponding 
conditional Liouville measure // on M can then be obtained as the contraction Lp{V) 
of the 2n-dimensional volume element V by a vector field {F(x) \ x G M} for which 
D F (H) = 1. Under our assumptions the phase space M is compact and, therefore, 
the measure \i is a finite, {^j-invariant Borel measure on M. The subject of this 
paper is the Hamiltonian flow with collisions (M, {ip 1 } , /j). 



We will frequently use another coordinate system (5h, 5v) in the tangent space 



For every interior point x = (q,p) of M we define the co dimension-one linear 
subspace T x of the tangent space T X ~M. of M at x as follows: 



M = {(q,p)eN\H(q,p) = l}. 



Factorization with respect to the flow direction 



T X N: 



(2.7) 




(2.8) 




(0,...,0; -1,...,-1) 



It is clear that the subspace T x of 7^ M is transversal to the velocity vector 

V(x) = U\x) 

(in (5h, 5v) coordinates) of the flow t G R}. 

Since the cu-orthocomplement of the tangent vector V(x) (x G intM) in X^N is 
precisely the space 7^M, we infer that the 2-form wfT^M naturally descends to 
the factor space 

T X M/ {XV(x)\ A G R} = T X M/V X , 

and it is a nondegenerate 2-form on that space. Moreover, the composition of 
the inclusion T x <— > T^M and the projection T^M — > T x ^AjV x provides a natural 
identification i : T x — > T^M/V^. The linearization (derivative) D^ f (x) : T X NL — > 
7^,t( x )M maps the line onto V^t^) and, therefore, it descends to a mapping 

: T X M/V X — > T^ (l) M/^, w . 

Thus, by using the abovementioned identification we can (and will) think of 
Di/; t (x) as a mapping D^i t {x) : T x — > %pt( x y The facts that the collision maps 
<E>i : <9M~ — > <9M+ are symplectomorphisms and {^*| t G R} is a Hamiltonian 
flow between collisions together imply that the linearization Dip t (x) : T x — > T^t^ 
preserves the nondegenerate 2-form uj. 

Characteristic (Lyapunov) exponents 

In the manifold N we introduce the Riemannian metric compatible with the 
2-form u> = J27=i dhi A dvf. 

n 

(2.9) ||(^,HI| 2 = E[(^) 2 + (^) 2 ]' 

i=i 

(5h, 5v) G T X N. It is worth noting here that there could be several other natural 
possibilities of defining a Riemannian metric in N, for example by using \\Sq\\ and 
\\Sp\\ 2 instead of \\Sh\\ 2 and ||5i>|| 2 . However, it is an easy consequence of (2.7) that 
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the ratio between these metrics would remain between two positive constants so 
that those constants would only depend on the masses of the particles and, being 
so, such a change in the Riemannian metric would not affect the characteristic 
exponents. It is obvious that the space T x (x G intM) is just the orthogonal 
complement (with respect to the Riemannian metric of (2.9)) of the line V x in the 
space T X ~WL. 

Thanks to the arising inner product in the space T x , we can now speak about 
the adjoint linear map 

(Dip^x))* : 7^t (a .) — > T x 

of the linearization Difj f (x) : T x — > T^t^ x y Both mappings are symplectic (and, 
therefore, invertible) with respect to the restriction of the 2-form to. 

We will be formulating Oseledets' Multiplicative Ergodic Theorem for a suitable 
Poincare section \1/ (a discrete time version) of the flow (M, {if)*}, y). Namely, we 
take ^> : <9M + — > <9M + the first return map to the boundary <9M + with respect to 
the flow {ip*}. The natural invariant measure /j, of \1/ on <9M+ can be obtained by 
taking fio(A) := lim^o -fJ,(A e ), A C <9M+, where 

A e = {x eM\3y e A3t, 0<t <e such that x = i/>*(y)} . 

It is a classically known result that the mapping (<9M + , \I>, /j, ) has (almost every- 
where) nonzero Lyapunov exponents if and only if the flow (M, {ip*}, A*) has the 
similar property. The reason is that the Lyapunov exponents of the map \& are just 
the corresponding characteristic exponents of the flow multiplied by the average 
(with respect to /j, ) return time t(x) (x G <9M + ) to the boundary <9M + , see, for 
instance, Lemma 2.2 of [W(1985)]. But the ergodic average t(x) is a /Uo-almost 
everywhere positive function. 

We introduce the shorthand C x = T x dM_ + for the points x G <9M + . The tangent 
space C x is a Euclidean space with the restriction of the inner product (2.9). The 
linearization D^(x) : C x — > and its adjoint (D^(x))* : C-%( x ) — > C x are now 
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symplectic isomorphisms between C x and C^( x y The integrability condition 

/ ||D*(x)|| dfio(x) < oo 

JdM+ 

of the Oseledets' theorem is an easy consequence of (10) from [W(1990)-I] and the 

following simple facts: 
d d 

(i) — (Sh) = -r(Sv) = (between collisions), 

(ii) the return time from <9M + to <9M + is bounded. 

The conditions of Oseledets' celebrated Multiplicative Ergodic Theorem, [0(1968)], 
[R(1979)], are fulfilled by the invertible measurable cocycle (5M + , // ): 

Multiplicative Ergodic Theorem. For \x-almost every phase point x G <9M + 
the following limit exists: 

lim \(DV n (x))* oWfi)]^ := A(x). 

n^oo 

The invertible limit mapping A(x) : C x — > C x is a symmetric, positive, symplectic 
linear transformation. The logarithms of the eigenvalues of A(x) 

-Xu-l(x) < -Xn-2(X) < ■ < -Ai(x) < < Ai(x) < • • • < Xn-l(x) 

are called the characteristic (Lyapunov) exponents of the invertible measurable co- 
cycle (aM+,*,M,/i ). 

The invariant cone field 

The symplectic linear space T X ~N is the direct sum Vl © Vi of the Lagrangian 
subspaces 

J Vi = {(Sh, Sv) G T X N\ Sv = 0} 
\V 2 = {(Sh, Sv) G T X N\ Sh = 0}. 
In Section 4 of [L-W(1995)j Liverani and Wojtkowski introduced the nondegenerate 

quadratic form Q in T^N associated with the decomposition T^N = V\ © V 2 as 

follows: 

n 

(2.10) Q ((Sh, Sv)) = {Sh; Sv) = Sh,S Vi . 

i=i 
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The corresponding positive cone (sector) C x C T X N between the Lagrangian sub- 
spaces Vi and V 2 is then defined as follows: 

(2.11) C x = {(8h, Sv) G T X N\ Q ((Sh, 6v)) > 0} . 

We will use the restriction of the quadratic form Q to the space T x and the inter- 
section C x fl T x , also (a bit sloppily) denoted by Q and C x . It is worth noting here 
that the Q-orthocomplement of the line V x in T X N is precisely the space 7^M and, 
therefore, the form Q descends to a nondegenerate quadratic form (also denoted by 
Q) on the factor space T X M_/V X = T x . 

Wojtkowski shows in Section 4 of [W(1990)-I] that if mi > m 2 > • • ■ > m n , then 
for every t > the linearized mapping Difj t (x) : T x — > T^t^ is Q-monotonic, i. e. 
for every tangent vector y G T x one has Q [(Dip t (x)) (y)] > Q(y) or, equivalently, 
the cone field C is invariant: 

(2.12) (Dip*(x)) (C x ) c C^ t(x) 
for every t > 0, x G intM, ^>*(a;) G intM. 

Definition 2.13. Let (cc) 1 1 > 0} 6e a nonsingular, positive orbit, x G intM. 
VFe say that the cone field C is eventually strictly invariant along the orbit 
{-0* (x) 1 1 > 0} iff there exits a number to > such that 

(Dip t0 (x)) (C x ) C int(C^ 0(x) ) , 

or, equivalently, Q [(Difj to (x)) (y)] > for every ^ t/ G C x . 

A major result (Theorem 5.1) of [W(1985)] is the following one: 

Theorem on Nonzero Characteristic Exponents. // the cone field C is even- 
tually strictly invariant along {^*(x)|t > 0} for ^-almost every x G M, then all 
characteristic exponents \i(x) of the cocycle (<9M + , /iq) are different from 

zero for hq- almost every x G <9M + . 

In Section 3 we will just check the conditions of this theorem for the falling ball 
system introduced before. 
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3. The Strict Invariance of the Cone Field 
(Proof of the Theorem) 

In this section we will be studying the strict cone invariance along a non-singular 
trajectory 

{^{x)\t ER} = {( qi {t), . . . , q n {t); v^t),..., v n (t)) \teR} 

of the Hamiltonian flow with collisions {•?/>* 1 1 G R} introduced in Section 1. We 
will always assume that mi > rri2 > > ■ ■ • > m n > and t = is not 
a moment of collision. Through references it was shown in the previous section 
that the quantities (5h;5v) = (6h\, . . . , 6h n ; Svi, . . . , 5v n ) (where we always as- 
sume Yl7=i = X^=i <^ Vi = 0) serve as suitable symplectic coordinates in the 
codimension-one subspace T x of the tangent space T X M of M at the phase point 
x = (9i(0), ••• , <Zn(0); i?i(0), . . . , v n (0)). Recall that the linear space T x is transversal 
to the flow direction and the restriction of the canonical symplectic form 

n n 

(3.1) uj = dhi A dv i = dqi A dpi 

i=l i=l 

of M is nondegenerate on T x . We also recall from the previous section that the 

1 p 2 
individual energy of the i-th particle is hi = rriiqi + -niivf = rriiqi + — L - and, 

therefore, dhi = midqi + rriiViSvi = rriiSqi + Vidpi. In his article [W(1990)-I] Wo- 

jtkowski introduced the nondegenerate quadratic form Q ((5h; Sv)) = (5h;5v) (the 

usual Euclidean inner product of 5h and dv) and the corresponding positive cone 

field 

(3.2) C x = {(Sh; 6v) E T x \ Q ((6h; 6v)) > 0} , 

see also the introduction of the present paper. He proved in the mentioned article 
that the cone field C x is invariant under the linearization of the flow {tfj t \ t e M}, 
i. e. Dipt^x) C C^t( x ) for every t > 0. (For a detailed discussion of the algebra 
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and geometry of an invariant cone field, see §4-6 of [L-W(1995)].) Furthermore, 
Liverani and Wojtkowski also proved in §4-6 of [L-W(1995)] that the eventually 
strict invariance of the cone field (i. e. D^iCx) C int (C^t( x )) for some t = t(x) > 
and for jU-almost every x G M) implies that all but one characteristic exponents of 
the flow {ip^te R} are nonzero! The only exceptional zero exponent corresponds 
to the flow direction. 

It follows from Wojtkowski's arguments between the Theorem and Proposition 1 
of Section 5 of [W(1990)-I] that in order to check the eventually strict invariance of 
the cone field along the studied nonsingular trajectory {^(x)] t G M} it is enough 
to prove that 

(A) for every vector ^ (0; Sv) G T x there exists a t > such that 

Q [£>V*((0;H)] > 

and 

(B) for every vector ^ (Sh; 0) G T x there exists a t > such that 

Q [Dip* ((6h;0)j\ > 0. 

Moreover, Wojtkowski's mentioned arguments from Section 5 of [W(1990)-I] ac- 
tually contain the proof of (A), provided that mi > mi > • • • > m n > 0. Here 
we briefly review his ideas. The formula (13) of [W(1990)-I] says that a tangent 
vector of the form (0; 8v) does not get changed at all by the linearization D<E> of 
the collision map $o corresponding to the floor collision. Suppose that a collision 
of type (i, i + 1) occurs at time tk- Denote the corresponding collision map by $j 
as in [W(1990)-I]. Suppose for a while that mi > mj+i. The equations (10) and 
(11) of the above article say that after pushing the tangent vector (0; Sv) through 
the collision + 1), the value of the Q form on the image D&i ((0;5v)) either 
becomes positive, or 5v~(tk) = Sv~ +1 (tk) and 5v + (tk) = Sv~(tk)- 

On the other hand, it also follows from (10) of the mentioned paper that in the 
case rrii = rrii+\ the linearization D$j of the collision map $^ simply interchanges 
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5vi and 5v i+1 : 5vf(t k ) = Sv~ +1 (t k ), 8vf +1 (t k ) = Sv~(t k ), 6vf(t k ) = Sv~(t k ) for 
j 7^ i,i + 1, and Sh + = Sh~ = 0. Therefore, it is quite reasonable to re-label 
the particles dynamically at every such collision by simply interchanging the labels 
i and i + 1. This is equivalent to allowing the particles with equal masses to 
freely penetrate through each other precisely the same way as Wojtkowski did in 
[W(1990)-I]. Then, as long as the Q form remains zero on the images of (0; 5v) : the 
images of (0; 5v ) under the linearization of the flow G M} remain the same, 

and dvi = Svj if the particles i and j collide on the considered trajectory segment. 
Since each particle i with rrii < mi must eventually bounce back from a heavier 
particle, and 1 is the sole heaviest particle by our assumption mi > m2, we obtain 
that every dvi must be the same. By our convention J27=i = ® ( we are a l wa Y s 
dealing with vectors from T x ), however, we infer that 5v = 0, provided that all 
future images of the considered tangent vector (0; Sv) have zero Q form. 

Thus, in order to prove the Theorem, it is enough to show that (B) holds true 
for (U-almost every x G M. This is what we are going to do. 

We begin with the definition of the "neutral space" M x of the nonsingular phase 
point x. (To be more accurate, Af x is going to be the neutral space of the positive 
orbit {ip t (x)\t > 0}.) The linear subspace M x of T x will be the precise analogue 
of the neutral space Mq (S^^x) of a positive orbit in a semi-dispersing billiard, 
originally and essentially introduced by Chernov and Sinai in [S-Ch(1987)], and 
later heavily used by Kramli, Szasz and myself in the several proofs of ergodicity 
for hard ball systems, [K-S-Sz(1991)], [K-S-Sz(1992)], [Sim(1992)-I-II], [S-Sz(1995)], 
[S-Sz(1996)]. 

Definition 3.3. 

M x := {(Sh;0) G T X \Q [Djj* ((5h;0))] =0 Vt > 0} . 

It is easy to convince ourselves that, indeed, Af x is a linear subspace of T x . The 
main result of this section (immediately proving the theorem) is 
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Main Lemma 3.4. For fi-almost every (nonsingular) phase point x G M we have 

Proof. The proof will be based on a few lemmas. Denote by < t\ < t 2 < . . . , 
tk /" oo the sequence of all collision moments on the positive orbit t > 0}. 

(These collision moments do not have a finite accumulation point, see the Appen- 
dix.) It is obvious that the tangent vector Dip 1 ((5h;0)) := (Sh(t); Sv(t)) does not 
change between collisions (after the natural identification of the tangent spaces of 
M at different points). Between the Theorem and Proposition 1 of Section 5 of 
[W(1990)-I] Wojtkowski proved that for any neutral vector (Sh; 0) G M x we neces- 
sarily have that 

(1) 8v{t) = 0, i. e. ((5h; 0)) = (5h(t); 0) W > 0; 

(2) 5h~(tk) = Sh + (tk) and 8h\{tk) = if the first particle bounces back from 
the floor at time tk-, i. e. q\{tk) = 0; 

(3) 5h + (tk) = R*Sh~(tk) if tk is a moment of an (i, collision (1 < i < n — 1), 
where Ri is the following n x n matrix: 



(3.5) 



Ri 



I 

o ■•• 

... 

... 







ll 1 - li 

1 + It -H 







... 

... 

'•• 

1 



Here 7^ = — — 7711+1 anc j ^j^g f our entries containing 7, fill up the intersections of 

mi + m i+1 

the i-th and j-th rows and columns. (For property (3), see also (10) of [W(1990)-I].) 

Lemma 3.6. For every vector (Sh; 0) G Af x we have J27=i 3hiVi(0) = and, hence, 
T l 7=iSh i (t)v i (t) = 0for allt>0. 

Proof of Lemma 3.6. Consider the quantity 

n 

w(t) :=^qi(t)6hi(t) 
i=i 
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well defined for all t > 0. The obvious relation ^-5hi(t) = implies that between 
collisions 

j n 

(3.7) -wit) = Vi(t)Shi(t) = (v(t); Sh(t)} 



dt 

i=i 



and 



(3.8) ^w(t) = -J2^i(t) = 0. 

i=l 

(Here we took advantage of the fact that in = — 1 is the gravitational acceleration.) 
Thus w(t) is a linear function of t between collisions. It is a straightforward con- 
sequence of (2) that the function w is continuous and it even does not change its 
slope at a collision with the floor. If, on the other hand, a collision of type (i, i + 1) 
takes place at time tk (1 < i < n — 1), then we have that = qi+i(t k ) and 

the compound velocity vector v(t) gets transformed by the matrix Ri at time tk' 
v+(t k ) = RiV-(t k ), see (9) of [W(1990)-I]. These facts and property (3) imply that 
w + (t k ) = w~(t k ) and 

jw{t k + 0) = (Sh(t k + 0); v(t k + 0)) = (R*5h(t k - 0); R iV {t k - 0)) = 

(5h(t k - 0); Rjv(t k - 0)) = jw{t k - 0). 

Here we used the obvious relation Rf = 1. We have seen therefore that w(t) is an 
(inhomogeneous) linear function of t > 0. 

Sublemma 3.9. TTie function w(t) is bounded (t > 0). 

Proof. It is enough to prove that the quantity Yl7=i ~~ t^^i(^)] 2 i s a constant 

Tfli 

function of t, since the positions qi(t) are obviously bounded. However, the quantity 

n 1 n 

^ — [^(t)] 2 = ^^(t)%(t) = (6h(t);8q(t)) 



mi 

%=i % i=i 



is constant between collisions and, by (2), it does not change its value at a floor 
collision. Furthermore, if a collision of type (z,z + 1) takes place at t k , then the 



16 

vector 5q(t) = D~ 1 5h(t) - where D = diag(mi, . . . , m n ) is the diagonal matrix 
with the masses as entries - gets transformed by the matrix DR*D~ X = Ri, i. e. 
5q(tk + 0) = Ri5q(tk — 0). Therefore, according to (3), we see that 

(Sh(t k + 0);Sq(t k + 0)> = (R*6h{t k - 0); RM{t k - 0)> = (Sh(t k - 0);Sq(t k - 0)> . 

(The above arguments are precisely the arithmetic background of the conservation 
of the kinetic energy at an (i,i + 1) collision.) 
Hence the sublemma follows. □ 

The assertion of the sublemma, together with (3.7), now proves Lemma 3.6. □ 

Finishing the Proof of Main Lemma 3.4 

Set 

(3.10) Xa = {i£ M| x is nonsingular & dimA/^ = a 7 } 

for d = 0, 1, ... ,n — 1. We want to show that for every d > the set has 
(i measure zero. Fix a number d > and an arbitrary phase point xq G X^- 
We will prove that x has a suitably small, open neighborhood U in M with 

n(x d nu ) = o. 

First choose a number r > with the following properties: 

(i) t is not a moment of collision in the positive trajecory of xq; 

(ii) M X0 = M({^ t (x )\0 < t<r}), i. e. the relation (Vt, < t < r) 
Q [Difjt^Sh; 0))] = implies that (Vt > 0) Q [D^{{8h; 0))] = 0. 

It is a very important consequence of properties (l)-(3) above that the neutral 
space of a finite orbit segment is completely determined by the symbolic collision 
sequence of that segment, i. e. by the types of collisions! Therefore, one can surely 
find a small, open neighborhood Uq of xq in M such that for every element y G Uq 
we have that 
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(i)' the symbolic collision sequence of {'0 t (y)| < t < t} is the same as of 

{V>*M|0<t<r} 

and r is not a collision moment of these orbit segments; 

(nyX({^(y)\0<t<r})=X Xo . 
Then we conclude that 

(dimM X0 =d (>1), 

I (VyeC/o) M y <zM XQ . 
Therefore, we get that M y = M XQ for every y £ Uo (~) Xa- According to Lemma 

3.6, however, the compound velocity v y of such a point y = (q y ,v y ) G Uq fl Xd 

is necessarily perpendicular to the Sh part of every neutral vector (8h; 0) G A/"a; . 

Since A/" Xo 7^ {0} and the velocities t> y of all points y £ U fill out an open subset of 

R n , we see that, indeed, fx (Uo fl X^) = 0. 

This finishes the proof of Main Lemma 3.4 and, therefore, the proof of the 

Theorem, as well. □ 

APPENDIX 
Degenerate Orbits 



We begin with a simple observation. Assume that we are given a phase point 
x = (q,p) G M with the property that there are k particles (1 < k < n) stuck on 
the floor with zero energy, i. e. q± = • • • = q^ = 0, p± = • • • = Pk = 0. Then the 
only natural way of defining the collisions of the orbit {^(x)] t G M} is such that 
these particles will remain standing still forever on the floor with zero energy. We 
call these trajectories the degenerate ones. 

Proposition A.l. Suppose that the trajectory {ifj 1 (x)\ t E R} is nondegenerate. 
Then the collision moments of this orbit can not accumulate at any real number t. 
In other words, there can only be finitely many collisions in finite time. 
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Proof. Assume the opposite, i. e. that there is a (say, positive) accumulation point 
of collision moments. Denote by to the smallest one of such positive accumulation 
points. Then all collision moments of the open interval (0, to) form an increasing 
sequence < t± < t<i < ■ ■ ■ < to such that lim m ^oo t m = to- Denote, as usual, 
by o~k = (ik,ik + 1) (0 < ik < n — 1) the type of the collision taking place at 
time tk, where (0, 1) means the collision with the floor. We say that an integer i 
(0 < i < n — 1) is essential if and only if there are infinitely many natural numbers k 
with i k = i. By classically known results ([G(1978)], [V(1979)], [G(1981)]), without 
the floor collision there can only be finitely many collisions in finite time. Therefore, 
the set of essential indexes i has the form {0, 1, . . . , a}, < a < n — 1. We can 
assume that the origin t = is already chosen in such a way that ik < a for every 
positive integer k. We are now focusing on the limit \im t y to Vi(t) of the velocities 
Vi(t) of the particles as t / to. 

Lemma A. 2. For every i, 1 < % < a + 1, the limit 

lim Vi(t) := v-(t ) 

exists. 

Proof. We start with the case i = a + 1. Since v a +i = — 1 between collisions and 

v a+i(tk) > v~ +1 (tk) for Zfc = a, we conclude that 

oo 

(A3) lim v a+1 (t) = v a+1 (0) -t + J2 K+i(tk) - v~ +1 (t k )] . 

t/t0 k=i 

This settles the case i = a+1. 

Suppose now that 1 < i < a and the lemma has been proved for i + 1, . . . , a + 1. 

Then again we have Vi = — 1 between collisions of type i — 1 or i and vf(tk) > v~ (t^) 

if ik = i — 1, while vf(tk) < v~ (tk) if ik = i- Since the lemma is supposed to be 

valid for i + 1, . . . , a + 1, one concludes that 

[ v r(tk)-Vi(tk)] <oo. 

&i = i 
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Thus, an argument similar to the one yielding (A. 3) provides 




Hence the lemma follows. □ 

Finally, since there are infinitely many collisions of type 1) (0 < % < a), 

we conclude that v~(to) < v~ +1 (to) (following from the fact that vf(tk) < vf +1 (tk) 
for ik = i) and, similarly, v~ +1 (to) < v~(to). (Here we use the natural convention 
vo(t) = 0.) Thus v~ (to) = for i = 1, . . . , a + 1 and, therefore, limt y to qi(t) = 0, 
as well. 

However, the subsystem {1, 2, . . . , a + 1} can not just loose its positive energy! 
The obtained contradiction finishes the proof of Proposition A.l. □ 

Remark. As it has been shown in [B-F-K(1995)], the number of collisions in a unit 
time interval is bounded in any semi-dispersing billiard. Since a small perturbation 
of a degenerate trajectory clearly provides an arbitrarily high number of collisions 
in unit time, we conclude that there is no way to introduce any Riemannian metric 
in the configuration space Q in such a way that 

(i) all sectional curvatures of that metric are nonpositive, 

(ii) the smooth components of the boundary <9Q are convex from inside Q, and 

(iii) the flow {V'*} is induced by the geodesic flow in Q and by the usual law of 
reflections at the boundary <9Q. 
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